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Abstract 

Semileptonic heavy — ► heavy and heavy — > light meson transitions are studied 
as a phenomenological application of a heavy-quark limit of Dyson-Schwinger 
equations. Employing two parameters: E, the difference between the mass 
of the heavy meson and the effective-mass of the heavy quark; and A, the 
width of the heavy- meson Bethe-Salpeter amplitude, we calculate for 
all decays on their entire kinematically accessible t-domain. Our study favours 
f B in the range 0.135-0.17 GeV and with E = 0.44 GeV and 1/A = 0.14 fm 
we obtain fS n (0) = 0.46. As a result of neglecting l/m c -corrections, we 
estimate that our calculated values of p 2 = 0.87 and f£ K (0) = 0.62 are too 
low by approximately 15%. However, the bulk of these corrections should 
cancel in our calculated values of Br(D — ► nlv) / Br{D — > Klv) = 0.13 and 

f5*(p)/f2 K (p) = 1.16. 
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I. INTRODUCTION 



Semileptonic meson decays are simple, experimentally accessible and only have a single 
hadron in the initial and final state. They are flavour- changing weak interaction processes 
and hence can be used as a means of extracting the elements of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, which characterise the difference between the mass eigenstates 
and the weak eigenstates in the Standard Model. For example, the K + — > 7r°e + z/ e (iQj) 
and K° L — > 7r ± e T z/ e (K® 3 ) decays currently provide the most accurate determination of 
|V^ S |(= 0.2196 ± 0.0023) JTJ, which would have been sin# c in the Cabibbo theory of weak 
interactions. The mechanism of the weak interaction is well understood. Hence, like elastic, 
electromagnetic form factors, these decays can also be used as a tool to probe the structure 
of the hadrons in the initial and final state. 

A major goal of current 5-meson experiments is to determine accurately the matrix 
elements V c b and V u b, the first of which is accessible via B — > D(D*)iv decays and the 
latter via B —>■ it(p)£v. The decays with a pseudoscalar meson in both the initial and final 
state are the simplest to study theoretically because they are only sensitive to the vector 
coupling of the W^-boson to the quarks and only two form factors are needed for a complete 
description. However, experimentally those with a vector meson in the final state provide 
the best statistics because the decay can proceed via both S- and D-waves. 

The B — > D(D*)£u decays proceed via aft^c transition and experimentally this is the 
closest one can come to realising a "heavy — > heavy" transition [T2J. It is in the analysis of 
these decays that heavy-quark symmetry [||]; i.e., an expansion of observables in Aqcd/"V, 
where m/ is the current-quark mass of the f — b,c quark, is most likely to be of use. However, 
in reality the AQ CD /m c -corrections, in particular, may nevertheless be large (~ 30)% and 
difficult to estimate in this case. 

Heavy Quark Effective Theory (HQET) || provides a systematic method for explor- 
ing the consequences of heavy-quark symmetry. It can be used to reduce the number of 
independent form factors required to describe heavy — > heavy decays, relating them to a 
minimal number of so-called "universal" form factors. However, it can't be used to calculate 
the g 2 -dependence of the form factors. This depends on the internal structure of the heavy 
mesons and its calculation requires the application of nonperturbative techniques. One such 
technique and its application to the calculation of these form factors is our focus herein. 

The methods of HQET are also not directly applicable to the decays B — > 7i(p)£u£, 
D — > Kivi and D — > n£i>£, all of which have light mesons in the final state. A primary 
impediment is that the current-quark mass of the s-quark m s ~ 0(Aq CD ), hence AQ CD /m s 
is not a suitable expansion parameter. In addition, a theoretical description of these decays 
requires a good understanding of light quark propagation characteristics and the internal 
structure of light mesons. This is provided by the extensive body of Dyson-Schwinger 
equation (DSE) studies |,§ in QCD. 

The DSEs are a system of coupled integral equations whose solutions, the n-point 
Schwinger functions, are the fully-dressed Euclidean propagators and vertices for the theory. 
Once all the Schwinger functions are known then the theory is completely specified. To 
arrive at a tractable problem one must truncate the system at a given level. Truncations 
that preserve the global symmetries of a field theory are easy to implement 0. Preserving 
the gauge symmetry is more difficult but progress is being made 0. 
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In a general covariant gauge the dressed-gluon 2-point Schwinger function (Euclidean 
propagator), D fiu (k), is characterised by a single scalar function, which we denote Q{k 2 )/k 2 . 
Important here is the particular, qualitatively robust result of studies of the DSE for D^ v {k) 
that Q(k 2 )/k 2 is strongly enhanced in the infrared; i.e, its behaviour in the vicinity of k 2 = 
can be represented as a distribution |||9|]. The infrared enhancement in D^ik) becomes 
prominent for k 2 ~ 1 Ge V 2 and is not peculiar to covariant gauges [|10 ] . 



The dressed-quark propagator for a quark of flavour / can be written in the general 
form[] 

where Zf(p 2 ) is the momentum-dependent wavefunction renormalisation and Mf(p 2 ) is the 
momentum-dependent quark mass-function. The dressed-gluon propagator is an important 
element in the kernel of the DSE satisfied by Sf(p). In existing studies of this DSE that 
employ a dressed-quark-gluon vertex that is free of kinematic, light-cone singularities, the 
infrared enhancement in D lll ,{k) is sufficient to ensure that S(p) doesn't have a Lehmann 
representation. This entails the absence of coloured quark states from the spectrum; i.e., 
quark confinement [11]. If Q(k 2 ) < oo at k 2 = it is possible to obtain a solution, Sf(p), of 



the quark DSE that has a Lehmann representation |T2 . 

There is another important consequence of the infrared enhancement in Q(k 2 )/k 2 . The 
enhancement is characterised by a mass-scale u ~ Aqcd and for light quarks; i.e, u-, d- and 
s-quarks for which rrif < Aqcd, it generates a significant enhancement in Mf(p 2 ). A single, 
indicative and quantitative measure of this enhancement in Mf(p 2 ) is the ratio Mf/rrif, 
where Mf is the Euclidean constituent-quark mass defined as the solution of p 2 = M 2 (j> 2 ).[] 
The results 



*&~150, ^-10 (2) 



m u4 

demonstrate that the infrared enhancement in Q{k 2 )/k 2 leads to at least an order-of- 
magnitude infrared enhancement in Mf(p 2 ). It is nonperturbative in originP] and has impor- 
tant qualitative and quantitative implications for light meson observables, as illustrated in 
Refs. 1,0. 



The effect of the infrared enhancement in Q (k ) / k on M c ^{p ) is much less dramatic [14 



—^-1.5, ^-2.0. 3 

m b m c 



1 We employ a Euclidean metric with 5^ = diag(l, 1, 1, 1), 7^ = 7 M and {7^,71,} = 25^ u . A 
spacelike 4- vector, has k 2 > 0. 



2 Quark confinement entails that there is no "pole-mass" [11|, which would be the solution of 
p 2 + M 2 (p 2 ) = 0. 

3 The renormalisation-point-dependence of the current-quark mass affects the actual value of the 
ratio Mf /mf but not the qualitative features of this discussion. 
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In this case rrif ^> Aqcd and the momentum-dependence of M Cjb (p 2 ) is primarily perturbative 



in origin. As observed in Ref. |TJ| it is therefore a good approximation to write 

M b (p 2 ) = const := M b « M b E , (4) 

for p 2 ^> — m 2 , although the 6-quark is still confined and there is no pole mass. For the 
same reason Z b (p 2 ) = 1 is also a good approximation. This and Eq. ([|) form the basis of 
the heavy-quark limit of the DSEs explored in Ref. [TJj wherein, on the domain explored by 



heavy — > heavy semileptonic decays, the dressed-6-quark propagator was approximated by 

S b (p) = , 1 • (5) 
ij ■ p + M b 

In Ref. |D| the dressed-c-quark propagator was approximated by an analogous expression: 

S c (p) = , 1 • (6) 
Z7 • p + M c 

However, the justification of this is less certain because the momentum dependence of Z c (p 2 ) 
and M c (p 2 ) is significantly more rapid. The approach employed in Ref. |TJ|] is one means of 



is 



exploring the fidelity of this approximation, as are the direct studies for which Ref. |14 
the pilot. 

Our aim herein is a unified description and correlation of semileptonic heavy — > heavy 
and heavy — > light meson transitions as an extension of the application of DSE methods. 



We follow Ref. |T5|] in describing the b- and c-quark propagators by Eqs. (^) and ([]), respec- 
tively, and in our analysis we consider the effects and limitations of Eq. (||). These equations 
represent the primary, exploratory hypothesis in our study because the propagation char- 
acteristics of light-quarks and the structure of light meson bound states is well understood 
following the extensive application of DSE methods in this domain |5|,|ll||l7j . In Sec. II we 
define our approximation to the matrix elements describing B(D) — > n(K)£i> transitions 
and fully specify a heavy-quark limit of our DSE application. Our results are presented and 
discussed in Sec. Ill and we make some concluding remarks in Sec. IV. 



II. SEMILEPTONIC DECAYS 



Our primary focus is the pseudoscalar — > pseudoscalar semileptonic decay 

PhM^PhMZv, (7) 

where Pjj t represents either a B or D meson with momentum pi (p 2 = —m 2 Hi ) and Ph 2 can 
be a D, K or 7r meson with momentum p 2 {p\ = —iti 2 h 2 )- The momentum transfer to the 
lepton pair is q := p\ —p 2 - A review of these decays is provided in Ref. [0] and a theoretical 
study of the light — > light transitions is presented in Ref. |T8|| . 
The invariant amplitude describing the decay is 

A(P Hl - P H M = % V qQ h,{\ - l5 )vM^ PH > ( Pl ,p 2 ) , (8) 
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where Gf is the Fermi weak-decay constant, V q Q is the appropriate element of the CKM 
matrix and the hadronic current is 

M^ PH >( Pl ,p 2 ) := (PhMIMIPhM) (9) 

= /+(*)(pi+M, + /-(*)<&,, (10) 

with t := — q 2 . The form factors, f±{t), contain all the information about strong interaction 
effects in these processes and their accurate estimation is essential to the extraction of V q Q 
from a measurement of a semileptonic decay rate: 



T(P Hl -> P H Ju) = J?L \v qQ \ 2 -±- jf~ <K |/+(t)| 2 [(*+ - *)(*- - 1)] 3/2 , (11) 



G 2 ? , „ , o 1 

with £± := (wijj-j ± vnn^) 2 and neglecting the lepton mass. 



A. Impulse Approximation 

In impulse approximation 

M^ PH >(p 1 ,p 2 )= (12) 

Mr r- n 

^ y d 4 fctr [r H2 (A;; -p 2 )S q (k + p 2 )iV^ Q (k + p 2 , fc + Pi)5 , Q (/c +Pi)r^ 1 (A:;p 1 )^(A;)] , 
where: r^/cjpi) is the Bethe-Salpeter amplitude for the Hi meson; 

f H2 (k; -p 2 y := C^T H2 (-k; -p 2 )C, C = 7 274, (13) 

and M* is the matrix transpose of M; and V^(fci,fc2) is the vector part of the dressed- 
quark-W-boson vertex. 

i. Quark Propagators 
The dressed quark propagators, Sf(p), in Eq. fll2"|) are the solution of 

S(p)- 1 = t 1 -p + m hm + J ^ g *D^{p-q)^S{q)T a v {q,p), (14) 



(27T) 

where: D^ u (k) is the dressed-gluon propagator; T"(g,p) is the dressed-quark-gluon vertex; 
mbm is the current-quark bare mass; and one can write Sf(p) in the general form 

s f (p) = -ij-p4(p 2 )+4(p 2 ), (is) 

which is completely equivalent to Eq. ([!]). A thorough discussion of the numerical solution 
of Eq. ([14]) , including a discussion of renormalisation, is given in Ref. [|17j . 
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a. Light quarks. Herein, for the light u-, d- and s-quark propagators, we do not directly 
employ a numerical solution of Eq. fll4) . Instead we use the algebraic parametrisations of 
these solutions developed in Ref. JT5[ because they efficiently characterise the essential and 
robust elements of the solution obtained in many studies 0] of the quark DSE: 

a f s (x) = 2m s T(2(x + rh 2 )) + J zr (6 1 x)J r (6 3 x) (b f + b{F{ex)) , (16) 



where: / = u,s (isospin symmetry is assumed), 

1 - e~ y 

Hv)--=—^—, (18) 

x = p 2 /(2D); fhf — rrif /\ / 2D; and 

4(x) := V2D a f s (p 2 ) , (19) 
4(x) :=2Da f v (p 2 ), (20) 

with D a mass scale. This algebraic form combines the effects of confinement and dynam- 
ical chiral symmetry breaking with free-particle (asymptotically-free) behaviour at large, 
spacelike p 2 . The parameters fhf, % 3 in Eqs. (|T6| ) and ( p!7[ ) take the values 

m f 4 b{ b{ b{ 

u : 0.00897 0.131 2.90 0.603 0.185 , (21) 
s : 0.224 0.105 2JH) 0.740 0.185 

which were determined in a least-squares fit to a range of light-hadron observables. The 
values of b\ 3 are underlined to indicate that the constraints b\ 3 = 6" 3 were imposed in the 
fitting. The scale parameter D = 0.160 GeV 2 . 

b. Heavy quarks. As described in Sec. |, and exploited in Ref. ||15|| , the momentum- 



dependence of Zf(p 2 ) and Mf(p 2 ) is much weaker for the heavy-quarks than it is for the 
light-quarks. This is illustrated for two different but related DSE-models in Refs. [14jl5H and 



justifies Eq. @ for the 6-quark and the cautious, exploratory use of Eq. @ for the c-quark. 

These equations provide the origin of heavy-quark symmetry in the DSE framework. Its 
elucidation is completed by introducing the heavy-meson velocity, v^, via 

Pi^:=m Hl v^:=(M fQ +E)v ll , (22) 

where v 2 = — 1 and E > is the difference between the heavy-meson mass and the effective- 
mass of the heavy-quark, M/„ . Equations and (|5|) then yield 




Exact heavy-quark symmetry arises from completely neglecting the l/Mf Q corrections in all 
applications. The mass of the 6-quark may justify this as a quantitatively reliable approx- 
imation but in making the same truncation for the c-quark one may expect quantitatively 
important corrections. 
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2. Bethe-Salpeter amplitude 



As discussed in Refs. |M],[L7||, the meson Bethe-Salpeter amplitudes in Eq. (|T2|) are the 



solution of the homogeneous Bethe-Salpeter equation: 

d 4 q 



lT H (k;P)} tu = J-^[XH(q;P)] sr Ki:( q ,k;P), (24) 
where 

XH (q; P) = S Q (q + P)T H (q; P)S q ,(q) , (25) 

Sf are the dressed-quark propagators, and r,. . . ,u represent colour-, Dirac- and flavour- 
matrix indices. In Eq. (^4j) K^(q, k; P) is the fully-amputated quark-antiquark scattering 
kernel. K^(q, k; P) is a 4-point Schwinger function obtained as the sum of a countable infin- 
ity of skeleton diagrams. It is two-particle-irreducible, with respect to the quark-antiquark 
pair of lines, and does not contain quark-antiquark to single gauge-boson annihilation dia- 
grams, such as would describe the leptonic decay of a pseudoscalar meson. The numerical 
studies of Ref. fl7j employed a ladder-like approximation: 



K™(q,k;P) = -g 2 D^(k- q ) U-) 7 „- , (26) 




which is consistent with the impulse approximation for M^ HiPh2 (pi,p 2 ) and is a quantita- 
tively reliable truncation for light, pseudoscalar mesons because of cancellations, order-by- 
order, between higher order diagrams in the skeleton expansion for K M. Ref. |bl is a first 



step in exploring the application of the methods of Ref. |L7j to mesons containing at least 
one heavy quark. 

a. Heavy meson Bethe-Salpeter amplitudes. Herein we do not use a numerical solution 
of Eq. (^4j) for the heavy-meson Bethe-Salpeter amplitude because we judge that our present 
studies are inadequate. One limitation, for example, is that simple ladder-like truncations 
do not yield the Dirac equation when the mass of one of the fermions becomes infinite and 
that defect may also be manifest in our study. Postponing the detailed exploration of this 
and other questions we employ instead an Ansatz motivated by the studies of Ref. |HJ and 
used efficaciously in Ref. p5|: 



T Hlf (k; Vl ) = 7 5 (l + ivy ■ v) V (k 2 ) , (27) 
where Mh 1} is the canonical Bethe-Salpeter normalisation constant. Using Eq. (|23| ) 

^, = — #3 rdu V (z) 2 (4(z) + yfr4{z)) :=—^, (28) 
f m H 1 f o2ti z jo v ' m Hlf Kf 

where z = u — 2Ey/u and / labels the light-quark flavour. 

In a solution of the Bethe-Salpeter equation the form of <p{k 2 ) is completely determined. 
However, here it characterises our Ansatz and as our primary form we choose 

(f(k 2 ) = exp (-k 2 /A 2 ) , (29) 
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where A is a free parameter. In studies of heavy — > heavy transitions Jl5| we found that, as 
long as f{k 2 ) is a non-negative, non-increasing, convex up function of k 2 , the results were 
insensitive to its detailed form. As we shall see below, through a comparison of the results 
obtained using Eq. (^) and those obtained with 



(p{k 2 



A 2 



(30) 



k 2 + A 2 

the same is true herein. Qualitatively, a primary requirement for an understanding of all 
the processes we consider is simply that the heavy meson be represented by a function that 
describes it as a finite-size, composite object: 1/A is a rough measure of that size. 

b. Light meson Bethe-Salpeter amplitudes. Just as for the light-quark DSE, there have 
been numerous studies @,[| of light mesons using Eq. (|24]) and a thorough discussion of the 
numerical solution, including a discussion of renormalisation, is presented in Ref. [O. The 
light, pseudoscalar meson Bethe-Salpeter amplitude has the general form 



T H (k;P) = l5 



iE H {k-P)+ 1 -PF H {k-P) 



(31) 



+ 1 -kk-P G H (k; P) + <v k^P v H H (k; P) 



Until recently it was assumed that in quantitative phenomenological applications one could 
neglect all but En(k; P) in describing the light, pseudoscalar meson and this was the assump- 
tion of Ref. 
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However, a systematic study of the quark DSE and meson Bethe-Salpeter 
equation |l7j demonstrates that the other functions are both qualitatively and quantita- 
tively important. A reanalysis of elastic form factors using all amplitudes and refitting the 
parameters characterising the quark propagators is therefore necessary. It is underway but 
incomplete |2T[ . 



Herein we use the parametrisation of the light meson Bethe-Salpeter amplitude deter- 
mined in Ref. [|19j and the results |I7j that it is a good approximation to treat En{k; P) = 
E H (k 2 ) and E 7T (k 2 ) = E K (k 2 ) := £{k 2 ); i.e., we use 

T„ = ^ K (k 2 ) =i lb £(k 2 ) 

y/2 Co e~ 



H=TT,K 



-k 2 /[2D] 



£{k 2 ) = ^ 



+ os(k 2 )\ 



M* a )|, 



(32) 
(33) 
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and therein yields the experi- 
is flawed. However, the light 



where the parameter Co = 0.214 GeV was fixed in Ref. 
mental value /„- = 0.131. For the kaon /# = 0.196 GeV. 

In principle, neglecting the other amplitudes in Eq. (| 
quark propagators of Eqs. fli~6|)-(|2lD were also fixed under this assumption and it is the 
combination of these parametrisations in Eq. ( p5|) that appears in the calculation of hadronic 
observables and reproduces available data. Therefore, if practiced judiciously, neglecting 
the other amplitudes can still provide quantitatively reliable results. To illustrate this we 
note that a preliminary reanalysis of the electromagnetic pion form factor |2l|], using all 
the amplitudes in Eq. (|3T| ) and refitting the w-quark propagator parameters in Eq. (pT|). 
yields results that are qualitatively indistinguishable from those obtained in Ref. [l!| for 
q 2 < 20 GeV 2 . It is only for q 2 > 20 GeV 2 that the qualitative and quantitative importance 
of F v and becomes manifest: these are the dominant amplitudes at large q 2 and ensure 
that q 2 F 7T (q 2 ) = const, up to lng 2 corrections. 



S 



3. Quark- W-bo son vertex 



Vp{k\,k2) in Eq. (|T^) satisfies a DSE that describes both the strong and electroweak 
dressing of the vector part of the quark- W-boson vertex. Solving this equation is a problem 
that can be addressed using the methods of Ref . [17] . However, we postpone this problem for 
the present and note instead that from this DSE one can derive a Ward-Takahashi identity 

(h - fc^iV^fo, k 2 ) = Stffa) - S h \k 2 ) - (m fl - m h ) Y^\h, k 2 ) , (34) 

where r{ 1 ^ 2 (fc 1 , k 2 ) is the scalar vertex, which in the absence of interactions is simply the 
diagonal unit matrix in Dirac space. This identity can be used to constrain the form of 
Vj^fe^ki, k 2 ), as the QED analogue has been used to constrain the dressed-quark-photon 
vertex [E2l. 



When fi and f 2 are both heavy quarks then the ability to neglect gluon dressing, as 
manifest in Eq. (|5D, entails 

(m fl - m h ) T^ih, k 2 ) « (M fl - M h ) 1 D . (35) 

This justifies the approximation, used efficaciously in Ref. [fT5| , 

V**(*i,A*) = 7„ (36) 

thereby amplifying the simplifications accruing in the heavy-quark limit. As demonstrated in 
Ref. Hl~8| , even in the case where both quarks are light, improvements to Eq. (|36| ) only become 
quantitatively significant (~ 10%) in the magnitude of f+(t) at the extreme kinematic limit: 
t = t-. Hence we use Eq. ( |36D in all calculations described herein. 



B. Semileptonic decays in the heavy-quark limit 



1. B f 



Di 



Using Eqs. flI5|), ( P3| ) and (P7D, we find JT5| from Eqs. ( |TUD and ( P^D that, at leading order 
in 1/rriH where uih is the heavy-meson mass, 



m D ± m B 

f±(t) = 2 ; H , 

y/m Df mB f 

^ ) = ^32^7o dT W L dULp{Zwf 



(37) 
(38) 



with W = 1 + 2t(1 - t)(w - 1), z w = u - 2EJu/W ancfl 



4 The minimum physical value of w is w m i n = 1, which corresponds to maximum momentum trans- 
fer with the final state meson at rest; the maximum value is w maiX ~ (m^ + m^^)/ (2mBfiriDf) = 
1.6, which corresponds to maximum recoil of the final state meson with the charged lepton at rest. 
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The canonical normalisation of the Bethe-Salpeter amplitude, Eq. (p8|), automatically en- 
sures that 

d f {w = 1) = 1. (40) 

Equation (|38D is an example of a general result that, in the heavy-quark limit, the semilep- 
tonic Hf — > H'j decays of heavy mesons are described by a single, universal function: 



2. Heavy —* Light 



Using Eqs. (0), fl23|) and (|27|) , and following the method outlined in the appendix, we 
find from Eqs. flit]) and (0) 



where 

F q ,(t;E,m Hl ,m H2 ) = 

with 

W q! { 1 ,u,u) = 2r 2 

♦(■■ 

1 

+ 



4 /-i 



9 32** 

c?7 



d 



a s^)-^- a v(^) - (T v {z x )—a q s {z 2 ) 



U V 



a u Jz 1 )aUz 2 )+uuaUz 1 )aUz 2 ) 



d 



+ (z x + u v Mfj 1 ) cry (zi) cr q v (z 2 ) - 2m 2 H2 r 2 a v (z x ) —a q v (z 2 ) 

az 2 



and 



z x = u 2 — 2uv E , 

z 2 = u 2 — 2u v (E — X) — m 2 H2 + 2i mn 2 7 u \J\ — v 2 
X=(m Hl /2) [l + (m 2 H2 -t)/m 2 Hi ], 

T = u \J\ — V 2 J 1 — 7 2 . 



(41) 



dv \ u 2 du(p(zi) £(zx)Wg'('j,i',u) , (42) 



(43) 



(44) 
(45) 
(46) 

(47) 

We note that because we have assumed isospin symmetry a u also represents a <i-quark and, 
to illustrate Eq. ([41]), the B° — > 7i~£ + V£ decay is characterised by 



f+ w (t) = «d^^2 F d {t- E, m B , m n ) 



(48) 
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C. Leptonic decays of a heavy meson 



We are also interested in the leptonic decay of a heavy, pseudoscalar meson, which is 
described by the matrix element 



(0\qi,l 5 Q\P Hl (p)} ■ = 
N c f 

fmPn = j^ti J d 4 ktr [n/5^S Q (k + p)T Hl (k;p)S q (k)} 



(49) 



where fu 1 is a single, dimensioned constant whose value describes all strong interaction 
contributions to this weak decay. For light mesons it has been studied extensively ||,r7| and 
with this normalisation f n = 0.131 GeV. Using Eqs. fli"5D , (p3|), (p7|) and (|2q) one obtains an 
expression for fu 1 valid in the heavy-quark limit : 



K 



Hi 



f 



y/m Hl 87T 2 



du {y/u - E) ip(z) a f s (z) + \yju cr v (z) 



(50) 



where z = u — lEy/u. It follows that in the heavy-quark limit 

f 1 

jHf OC 



\Jm Hf 



(51) 



This scaling law is counter to the trend observed in the light mesons, as highlighted 



in Ref. [14], where fn increases at least up to current-quark masses three-times that of 
the s-quark. Contemporary estimates of fjj and Jb, such as those analysed in Ref. 
suggest that Eq. (^1|) is also not obeyed by experimentally accessible heavy mesons. The 
determination of the current-quark mass at which the light meson trend is reversed, and 
that at which this heavy-quark scaling law is satisfied, is an interesting, open question. 



III. RESULTS AND DISCUSSION 

We have now defined all that is necessary for our calculation of the semileptonic heavy 
— > heavy and heavy — > light meson transition form factors and heavy-meson leptonic decay 
constants. We have two free parameters: the binding energy, E, introduced in Eq. ( |22| ) and 
the width, A, of the heavy meson Bethe-Salpeter amplitude, introduced in Eq. (|29|). The 
dressed light-quark propagators and light-meson Bethe-Salpeter amplitudes have been fixed 
completely in the application of this framework to the study of tt- and .fT-meson properties. 

Our primary goal is to determine whether, with these two parameters, a description and 
correlation of existing data is possible using the DSE framework. This was certainly true 



in our analysis of heavy — ► heavy transitions alone [ 15| . We found that the function £(w) 
necessarily has significant curvature and that a linear fit on 1 < w < 1.6 is inconsistent with 
our study. However, our calculated value of the slope parameter 



" 2 : = " 



(52) 

10 = 1 
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was too strongly influenced by the experimental fit to the B — > D data for that study to 
provide an independent prediction of p 2 .| Herein we eliminate this bias by excluding D- 
meson observables from our primary procedure for fitting E and A. This also facilitates an 
elucidation of where l/M c -corrections are important. 

Our key results are presented in column one of Table |. In obtaining these results we 
varied E and A in order to obtain a best, weighted least-squares fit to the three available 
lattice data points [B^j for f+ n and the experimental value |2B| for the B° — > -k~1 + v branch- 
ing ratio. In doing this we constrained our study to yield f B = 0.17 GeV from Eq. (p0|). 
which is the central value favoured in a recent analysis of lattice simulations |24{|, and used 



wb = 5.27 GeV. This fitting procedure assumes only that the 6-quark is in the heavy-quark 
domain; i.e., that l/M^-corrections to the formulae we have derived herein are negligible. 
Our calculated form of /+*"(£) is presented in Fig. [I]. A good interpolation of our result is 
provided by 



0.458 



/+(*) = 1 TT^ > m ™° = 5 - 67 GeV • ( 53 ) 



1 — t/m 



mon 



m r 



This value of m mon can be compared with that obtained in a fit to lattice data [25 
5.6 ±0.3. 

In Table H we compare our favoured, calculated value of f+ n (0) = 0.46 with this quantity 
obtained using a range of other theoretical tools. Since the t-dependence of /+*"(£) is an 
outcome of our calculation, the value we predict for f+ n (0) is the only one that allows 
simultaneous agreement between our calculations and existing results of lattice simulations 
and the measured branching ratio. If these data are correct then in our framework it is not 
possible to obtain a value of f+ n (0) that differs from this favoured value by more than 10% 
unless the calculated t-dependence is changed significantly. This could only be effected by 
a modification of the vertex Ansatz, Eq. flBlf), and hence the accuracy of our prediction can 
be seen as a test of the veracity of this Ansatz in the heavy-quark limit. 

In Fig. H] we present our calculated form of the function, £(u>), that characterises the 
semileptonic heavy — > heavy meson decays. We have compared our calculation with the 



experimental results of Ref . |28| and the following fits to the experimental data in Ref . |29[ : 

£(w) = 1 - p 2 {w - 1), p 2 = 0.91 ± 0.15 ± 0.16 , (54) 

, p 2 = 1.53 ±0.36 ±0.14. (55) 



= — 2 — exp 
w + 1 



(l-2p 2 )^ 
V H J w + 1 



Our calculated result for p 2 is close to that in Eq. (|5lD but our form of £(w) has significant 
curvature and deviates quickly from the linear fit. The curvature is, in fact, very well 
matched to that of the fit in Eq. (|55D, however, the value of p 2 listed in that case is very 
different to our calculated value. 

In Ref. ]T3] we fitted E and A to the nonlinear form in Eq. ( |55D and fitted it exactly. We 



believe that part of the discrepancy observed here is due to our neglect of l/M c -corrections 
in the calculation of £ (w) , the magnitude of which is exposed because of our newfound ability 



Til our framework the minimum possible value for p 2 is 1/3 JT 
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to constrain our parameters without referring to D-meson observables. Nevertheless, the 
agreement between this calculation and the data is reasonable, with the difference largest at 
w m£H where it is a little more than one standard deviation. Hence l/M c -corrections cannot 
be too large. 

In Fig. H we present our calculated form of f+ K (t). The t-dependence is well- 
approximated by a monopole fit. Our favoured, calculated value of f+ K (0) = 0.62 is 
approximately 15% less than the experimental value |TJ. We interpret this as a gauge of 
the size of l/M c -corrections. These corrections are expected to reduce the value of the 
D-meson leptonic decay constants from that obtained using Eq. Q5TD. A 15% reduction in 
the D-meson leptonic decay constants in column one of Table [I] yields fo = 0.24 GeV and 
fn s = 0.26 GeV, values which are consistent with lattice estimates |24j] and the latter with 
experiment [|30 . 

We have also calculated f+(t) and find that on the kinematically accessible domain, 



< t < (m D 



the following monopole form provides an excellent interpolation 



frit) 



0.716 



t/m 



2 

mon 



2.15 GeV. 



(56) 



We note that a naive vector meson dominance assumption would lead one to expect m Tl 
rriu* = 2.0 GeV . Using (E, A) from Table | we obtain 



Br(D 



nils 



Br(D -> Klv) 



2.47 



V, 



cd 



0.13 



(57) 



for \V c d/V cs \ 2 = 0.051 ± 0.002 [l|], and in this ratio the bulk of the l/M c -corrections should 
cancel. Experimentally 



Br(D c 



it e^v P 



0.111°;^ ±0.1 [1,H 



Br(D° 
Br(D + 



K~e + u P 



7T 



°e + u e ) 



Br(D- 



K°e+v P 



0.17 ±0.05 ±0.03 39 



(58) 
(59) 



We observe that if one makes the assumption of single-pole, D* and D* vector meson dom- 
inance for the t-dependence of the form factors and f+ K , respectively, one obtains the 
simple formula 



R n = 1.97 



/r(o) 



v cd 2 



(60) 



This approach has been employed |l[] in order to estimate f®'*^)/ f+ K (Q) = l-0jl o J ± 0.04 
or 1.3 ± 0.2 ± 0.1 from Eqs. @ and ©. We calculate 



f? K (o) 



1.16 



(61) 



It is incumbent upon us now to stress that we explicitly do not assume vector meson 
dominance. Our calculated results reflect only the importance and influence of the dressed- 
quark and -gluon substructure of the heavy mesons. This substructure is manifest in the 
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dressed propagators and bound state amplitudes, which fully determine the value of every 
quantity calculated herein. Explicit vector meson contributions would appear as pole terms 
in Vj}* 2 {k\, k 2 ), which are excluded in our Ansatz, Eq. fl36|). That simple-pole Ansatze 
provide efficacious interpolations of our results on the accessible kinematic domain is not 
surprising, given that the form factor must rise slowly away from its value at t — and the 
heavy meson mass provides a dominant intrinsic scale, which is modified slightly by the scale 
in the light-quark propagators and meson bound state amplitudes. Similar observations are 
true in the calculation of the pion form factor, as discussed in detail in Sec. 7.1 of Ref. 
and Sec. 2.3.1 of Ref. gOJ. 

In column two of Table | we present the results obtained when E and A are varied in 
order to obtain a best, weighted least-squares fit to: the lattice data on f+ n ; the B° — > 
tt~1 + v branching ratio; and the experimental data on reported in Ref. P8|. The latter 



introduce D-meson properties into our fitting constraints but their effect on our calculations 
is not very significant. The tabulated quantity most affected is the B° — > it~£ + i> branching 
ratio but this increases by only 15% and remains acceptably close to the experimental value. 
The effect that this modified fitting procedure has on the transition form factors is also 
small, as illustrated by the comparisons in Figs. [l]-|3[ Not surprisingly, the largest effect is a 
uniform 5% increase in the magnitude of f+ K (t). 

In Table |TTT| we present the results obtained using the different functional form for the 
heavy- meson Bethe-Salpeter amplitude in Eq. (|30D . A direct comparison with the results in 
Table [I] indicates that our results are insensitive to such details and hence are robust. The 
binding energy, E, is unchanged and the width, A, is smaller, as expected since Eq. ([30D 
does not decrease as rapidly with k 2 as the form in Eq. (p9|). A quantitative statement of 
this is that 

dk 2 (e" fc2 / A2 ) 2 = ±A 2 , (62) 



o 



DC 



A 



2 



2 



dk 2 I — I = A 2 (63) 
o \k 2 + A 2 J K J 

and A = 0.92 GeV ~ A/ y/2 = 1.0 GeV is just that reduction necessary to provide the same 
integrated strength for both amplitudes. 

Tables |V| and [V] provide a further elucidation of the impact of possible systematic errors 



in our calculation. These results are obtained through a repetition of the calculations that 
yield Table | but with constrained to be 0.135 and 0.205 GeV, respectively, which are 



the outer limits estimated in an analysis of contemporary lattice simulations f24 |. In the 



direct application of the methods of Ref. [T7|] to heavy mesons the value of fs would be a 
prediction. Herein, since we do not calculate but instead fit the heavy-meson Bethe-Salpeter 
amplitude, fs acts as a constraint on the width, A, of the Bethe-Salpeter amplitude, as seen 
in a comparison of Tables [|, fV| and [V|. The binding energy, E, is then the only true free 
parameter and it varies over a range of no more than 8%. Comparing these tables, we see 
that our results are not very sensitive to the value of fs in the range we have explored; i.e., 
our results are robust. 

We judge that the best description of the available data is obtained with f B = 0.17 GeV, 
with a lower value, — > 0.135 GeV, more acceptable than a higher one. The value of 
E = 0.44 GeV that provides this best description can be compared with the value of E hind ~ 
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0.25-0.35 GeV obtained in a lattice NRQCD simulation g|. The value of A = 0.92 GeV 
indicates that the heavy meson occupies a spacetime volume only 15% of that occupied by 
the pion. 



IV. CONCLUSIONS 

Using the same phenomenological Dyson- Schwinger equation (DSE) framework employed 



in successful studies of light meson observables as diverse as n-n scattering [42] and diffractive 



electroproduction of vector mesons [13], we have analysed semileptonic heavy — > heavy and 
heavy — > light meson transition form factors. In this application we introduced and explored 
a heavy-quark limit of the DSEs based on the observation that the mass function of heavy 
quarks evolves slowly with momentum. 

With two parameters: E, the difference between the heavy-meson mass and the effective- 
mass of the heavy quark; and A, the width of the heavy meson Bethe-Salpeter amplitude, we 
obtained a uniformly good, robust description of all available B — > ix data with a prediction 
for f+ n {t) on the kinematically accessible t-domain. In analysing B — > D, D — > K and 
D — ► 7i transitions we estimated that l/m c -corrections to our heavy-quark limit contribute 
no more than 15%. A significant feature of our study is the correlation of heavy — ► heavy 
and heavy — > light transitions and their correlation with light meson observables, which are 
dominated by effects such as dynamical chiral symmetry breaking and confinement. 

This study can be extended, with the application of the framework to semileptonic decays 
with vector meson final states using no additional parameters. It can also be improved, for 
example, by an exploration of the effect of more sophisticated Ansatze for the dressed-quark- 
W-boson vertex and of the inclusion of all amplitudes in the light-meson Bethe-Salpeter 
amplitude with refitted light-quark propagators. 

A more significant qualitative improvement is the direct study of the Bethe-Salpeter 
equation for heavy mesons using the methods of Ref. |T7|]; Ref. |14| is the pilot. This pro- 
gramme involves the important step of critically analysing the reliability for heavy quarks 
of ladder-like truncations of the dressed- quark- ant iquark scattering kernel in both the quark 
DSE and meson Bethe-Salpeter equation. Addressing this question and developing an ef- 
ficacious truncation will allow a correlation of heavy- and light-meson observables via the 
few parameters that characterise the behaviour of the quark-quark interaction in the non- 
perturbative domain; i.e., relate both heavy- and light-meson observables to the long-range 
part of the quark-quark interaction. 
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APPENDIX: A DERIVATION 

A typical integral arising in the detailed analysis of Eq. has the form 

d A k 1 



To simplify it we introduce a Laplace transform for the functions Z(k 2 ) and a([k — p 2 ] 



dsZ(s)e- sk2 



-u[k-p 2 ] 2 



z(k 2 ) - 

poo 

a(\k — p 2 ] 2 ) = / dua(u)e~ 
Jo 

and a Gaussian representation of the heavy-quark propagator: 

k ■ v — E Jo 

Inserting these identities we obtain 

00 00 00 d^k 

J = J ds Z{s) J du &{u) J da J — j- exp{— sk 2 — a{k ■ v — E) — u[k + P2} 2 } 



00 00 00 

= J ds Z{s) J du a{u) J da exp {aE — up 2 , + {up 2 + \av ) 2 / {s + w) j 
000 
d 4 k ( r 

—j- exp I — (s + u) + («P2 + \av)/(s + u) 



(Al) 

] 2 ): 
(A2) 
(A3) 

(A4) 



(A5) 
(A6) 



Shifting variables: k — > k — (up 2 + \av)/(s + u) and subsequently a — > (s + tt)a, yields 



J = J dsZ(s) J dua(u) 



00 



(A7) 



y da (s + u) exp | — (s + u)(|a 2 — ai?) + «aJ ^ — p| j> ^ 



7T 



ds Z(s) I dua(u) Ida exp <j — (s + u)(^a 2 — aE) — uaX ——pi} (A8) 



-(s+u)k 2 



SU 



S + U 



where X := — v ■ p 2 , Eq. (EBf). Making use of the identities 



s + u ' 



exp 



SU 



s + u 



v 2 



I s + u 



71 



J dr exp j— st 2 — u (t + sfp^j j , 



y/s + u y/ir 



4- [°° dte-^ t2 



(A9) 
(A10) 



we obtain 
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2 f°° 

J = - da I dr (All) 



7T JO 

poo 



/oo , 
dt Z(a 2 - 2aE + r 2 + t 2 ) a (a 2 - 2aE + 2aX + (r + ^) 2 + t 2 
-oo ^ 

Introducing spherical polar coordinates 

a = uu, (A12) 
r = u Vl - v 2 i , (A13) 
t = u Vl - v 2 yfi - 7 2 , (A14) 

with u G [0, oo), v G [0, 1] and 7 G [—1,1], we arrive at 

J = - f 1 J 1 = Z 1 <2v /°° it 2 Z( Zl ) (t(^) , (A15) 

7T 7-1 vl — 7^ JO JO 

where, using p 2 , = —m 2 H2 , 

z x = u 2 - 2uvE, (A16) 



Zi = v 2 — 2uv (E-X) -m 2 H% +2im H2 uVl - v 2 ^1 - 7 2 . (A17) 

This is recognisably of the form in Eq. (^T|). 

Structures more complicated than Eq. (Al) arise in deriving the complete form of F q / } 
however, they can all be analysed and simplified using analogues of the method illustrated 
above. 
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FIGURES 



3.0 - 




t ( GeV ) 



FIG. 1. Our calculated form of f^(t): solid line - column 1, Table |L dashed line - column 2, Ta- 
ble |[ For comparison, the data are the results of a lattice simulation an d the light, short-dashed 
line is a vector dominance, monopole model: f+(t) = 0.46/(1 — t/m^»), mg« = 5.325 GeV. The 
light, dotted line is the phase space factor |/f^(0)| 2 [(t+ - t)(t- - t)f /2 /(Trm B ) 3 in Eq. flu]), which 
illustrates that the B — > ireu branching ratio is determined primarily by the small-g 2 behaviour of 
this form factor. 
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FIG. 2. A comparison of our calculated form of £(w) with recent experimental analyses. Our 
results: solid line - column 1, Table |; dot-dashed line - column 2, Table |. Experiment: data 
points - Ref. p^] ; short-dashed line - linear fit from Ref. p9jj , see our Eq. (|54|); long-dashed line 
- nonlinear fit from Ref. [29|, see our Eq. (55). The two light, dotted lines are this nonlinear fit 
evaluated with the extreme values of p 2 : upper line, p 2 = 1.17 and lower line, p 2 = 1.89. 
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FIG. 3. Our calculated form of f+ (q 2 ): solid line - column 1, Table |; dashed line - column 2, 
Table | For comparison the light, short-dashed line is a vector dominance, monopole model: 
f + (q 2 ) = 0.74/(1 - q 2 /m 2 D *), m D: = 2.11 GeV. 
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TABLES 





DATA/ESTIMATES 


f B = 0.170 GeV 


(E, A) (GeV) 




( n a /i o i /i no \ 
(^U.44z,1.4Uo J 


(^U.4oo,1.4Uo ) 


E 2 /N 




0.48 


1.22 


J + ^i.auev ) 


ns2 + ni7 ll?5ll 


0.84+ 


0.89* 


J + I, 1 ' v j 


J- • J. iJ _1_ \J • Zj O Li. fJII 


1.02 t 


1.09* 


f B7T (20 9 GeV 2 ! 


-L . O i7 —1— • (J KJ (Jll 


1.30 t 
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RrfR _> 77--/+;/") 


h gin4_|_rio_|_nol x i n-4 ira 


2.0 xl0 _4 + 


2.3xl0 _4 + 


J + l u J 


n 1 8 _> o 49 [ETll 


0.46 


0.48 




0.74 ± 0.03 Q 


0.62 


0.65 


efi net; i n (vit^ 


n oq _i_ n 1 o [?q| 

U.OO _l_ VJ.XO | | 


0.93 


0.93+ 


t(\\ i c _l n C\AK\ 


0.93 ±0.12 @ 


0.86 


0.86+ 


CM X.ZiO ZIZ U.UtJUl 


0.68 ± 0.06 || 


0.80 


0.79+ 


ffl ^fi + ft (WH 


0.66 ± 0.06 || 


0.76 


0.75 1 " 


pl-i Ar. _i_ n ()AK\ 

CM ±. L ±0 _l_ \J.\J'-tOJ 


0.58 ± 0.08 H 


0.71 


0.70+ 


p 


0.91 ±0.15 ±0.06 _ 
1.53 ±0.36 ±0.14 ^ J 


0.87 


0.92 


f Ba (GeV) 


0.195 ±0.035 [ 


n 


0.184 


0.184 


IbJIb 


1.14 ±0.08 |4 




1.083 


1.082 


fo (GeV) 


0.200 ± 0.030 [ 
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I 


0.285 


0.285 


/d s (GeV) 


0.220 ± 0.030 [ 


24 


I 


0.304 


0.304 


/d s //d 


1.10 ±0.06 |4 


1 


1.066 


1.066 



TABLE I. A comparison of our calculated results with available data when we require 
fs = 0.170 GeV, which is the central value estimated in Ref. and use Eq. (|29|). In each 

column the quantities marked by + are those used to constrain the parameters (E, A) by min- 



imising X 2 := J2iL\ ([Hi 3,10 ~ yf ata ]/cr(y)^ ata ) 2 , where N is the number of data items used. The 
results in the first column assume that heavy-quark symmetry is valid for the 6-quark but do not 
rely on this being true for the c-quark. We note that: 1) our values of fu and fjj 3 are obtained 
via Eq. (pl|) from fs and fs s , respectively, using ms = 5.27, ms s = 5.375, mn = 1.87 and 
rri£) s = 1.97 GeV; 2) the experimental determination of p 2 is sensitive to the form of the fitting 
function, e.g., see Ref. [29|; 3) an analysis of four experimental measurements of D s — > [iv decays 
yields f Ds = 0.241 ± 0.21 ± 0.30 GeV |C]. 
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Reference 


/f-(0) 


Our Result 


0.46 


Dispersion relations p7f 
Quark Model [31 1 
Quark Model [32| 
Quark Model f 33 1 

Light-Cone Sum Rules Q 

Quark Confinement Model [35 
Quark Confinement Model [36,37 


0.18 0.49 
0.33 ± 0.06 
0.21 ±0.02 
0.29 

f 0.29 direct 

1 0.44 pole dominance 

0.6 

0.53 



TABLE II. A comparison of our favoured, calculated result for fS n (0) with a representative 
but not exhaustive list of values obtained using other theoretical tools. More extensive and com- 
plementary lists are presented in Refs. [p7, 33,37]. 





DATA/ESTIMATES 


JB — U.l ( U Ucv 


(E,A) (GeV) 




(0.455,0.918) 


£ 2 /iV 




0.46 


/f "(14.9 GeV 2 ) 


0.82 ±0.17 


25j 


0.84t 


/f w (17.9GeV 2 ) 


1.19 ±0.28 


m 


1.02+ 


ff n (20.9 GeV 2 ) 


1.89 ±0.53 


S 


1.32t 


Br( J B° -» tt-^+i/) 


[1.8 ± 0.4 ± 0.3 ± 0.2] x 10" 4 |g] 


2.0 xl0" 4t 


/f*(0) 


0.18 -» 0.49 |3 


0.45 


/£*(<>) 


0.74 ± 0.03 § 


0.62 


£(1.085 ± 0.045) 


0.88 ±0.18 


n 


0.92 


£(1.18 ±0.045) 


0.93 ±0.12 


n 


0.84 


£(1.28 ±0.050) 


0.68 ± 0.06 


2§ 


0.77 


£(1.36 ±0.050) 


0.66 ± 0.06 




0.72 


£(1.45 ±0.045) 


0.58 ±0.08 


H 


0.67 




0.91 ±0.15 ±0.06 — 
1.53 ±0.36 ±0.14 ^ 


1.03 


f Ba (GeV) 


0.195 ±0.035 H 


0.180 


fBs/fB 


1.14 ±0.08 |24] 


1.061 


Id (GeV) 


0.200 ± 0.030 H 


0.285 


/d. (GeV) 


0.220 ±0.030 §§ 


0.298 


fc-Jfc- 


1.10 ±0.06 |24] 


1.044 



TABLE III. A comparison of our calculated results with available data when we require 
fs = 0.170 GeV, which is the central value estimated in Ref. [24], and use Eq. (|30|) . (See Table. | 
for additional remarks and an explanation of the symbols.) 



24 





DATA/ESTIMATES 


Sb = 0.135 GeV 


(E, A) (GeV) 




(0.457,1.138) 


/ r\ ACT' 1 1 or\ 

(0.466,1.135) 


/ 




0.50 


0.97 




n S9 + n 1 7 ifJ^ll 


0.86+ 


0.88+ 


f B7I 7l 7 Q HpV 2n I 


1 1 q x n oc 119 q| 


1.05+ 


1.08+ 


f Bn (90 Q fW 2N l 


1 on _i_ n c.o 119 ell 


1.36+ 


1.40+ 


R r f R° _^ -TT-f+lA 


[1 Q_|_n4_|_nQ_|_r)9l v 1 n— 4 loai 


2.1 xl0~ 4 + 


2.2xl0" 4 + 


J + l u J 


n is _> 49 [271 

• -L U ' V7 . t: i_/ lull 


0.46 


0.47 


J+ w 


n 74 _l n no 111 


0.64 


0.65 


£(1.085 ± 0.045) 


0.88 ±0.18 @ 


0.92 


0.92+ 




n no 1 n 1 9 119 oil 


0.85 


0.85+ 


p(-\ 90 1 n ncpA 


0.68 ± 0.06 H 


0.78 


0.78+ 


n 1.00 zt u.uouj 


0.66 ± 0.06 PS 


0.74 


0.73+ 




0.58 ±0.08 H 


0.69 


0.69+ 




0.91 ±0.15 ±0.06 — 
1.53 ±0.36 ±0.14 ^ 


0.96 


0.98 


f Bs (GeV) 


0.195 ±0.035 [ 




0.148 


0.148 


IbJIb 


1.14 ±0.08 [24 


1 


1.096 


1.096 


Id (GeV) 


0.200 ± 0.030 [ 


24| 


0.227 


0.227 


/d. (GeV) 


0.220 ± 0.030 [ 


24 1 


0.244 


0.244 


/d s //d 


1.10 ±0.06 §4 


I 


1.079 


1.078 



TABLE IV. A comparison of our calculated results with available data when we require 
Jb = 0.135 GeV, which is the lower bound estimated in Ref. pi]], and use Eq. (pS|). (See Ta- 
ble. I for additional remarks and an explanation of the symbols.) 
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DATA/ESTIMATES 


Sb = 0.205 GeV 


(E, A) (GeV) 




(0.469,1.677) 


(0.479,1.678) 


/ 




0.83 


1.45 




n S9 + n 1 7 ifJ^ll 


0.91+ 


0.94+ 


f B7I 7l 7 Q HpV 2n I 


1 1 q x n oc 119 q| 


1.11+ 


1.15+ 


f Bn (90 Q fW 2N l 


1 on _i_ n c.o 119 ell 


1.43+ 


1.49+ 


R r f R° _^ -TT-f+lA 


[1 Q_|_n4_|_nQ_|_r)9l v 1 n— 4 loai 


2.4 xl0~ 4t 


2.5xl0" 4 + 


J + l u J 


n is _> 49 [271 

• J. U ' V7 . t: i_/ lull 


0.49 


0.50 


J+ w 


n 74 _l n no 111 


0.66 


0.68 


£(1.085 ± 0.045) 


0.88 ±0.18 @ 


0.93 


0.93+ 




n no 1 n 1 9 119 oil 


0.86 


0.86+ 


p(-\ 90 1 n ncpA 


0.68 ± 0.06 H 


0.80 


0.79+ 


n 1.00 zt u.uouj 


0.66 ± 0.06 PS 


0.75 


0.75+ 




0.58 ±0.08 H 


0.71 


0.70+ 




0.91 ±0.15 ±0.06 — 
1.53 ±0.36 ±0.14 ^ 


0.89 


0.91 


f Bs (GeV) 


0.195 ±0.035 [ 




0.220 


0.220 


IbJIb 


1.14 ±0.08 [24 


1 


1.071 


1.071 


Id (GeV) 


0.200 ± 0.030 [ 


24| 


0.344 


0.344 


/d. (GeV) 


0.220 ± 0.030 [ 


24 1 


0.363 


0.363 


/d s //d 


1.10 ±0.06 |i 


1 


1.054 


1.054 



TABLE V. A comparison of our calculated results with available data when we require 
Jb = 0.205 GeV, which is the upper bound estimated in Ref. p4j , and use Eq. (p9|). (See Ta- 
ble. I for additional remarks and an explanation of the symbols.) 
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